Molecular dynamics in multidimensional space explains how mutations affect the association path of neomycin to a riboswitch.

Riboswitches are messenger RNA (mRNA) fragments binding specific small molecules to regulate gene expression. A synthetic N1 riboswitch, inserted into yeast mRNA controls the translation of a reporter gene in response to neomycin. However, its regulatory activity is sensitive to single-point RNA mutations, even those distant from the neomycin binding site. While the association paths of neomycin to N1 and its variants remain unknown, recent fluorescence kinetic experiments indicate a two-step process driven by conformational selection. This raises the question of which step is affected by mutations. To address this, we performed all-atom two-dimensional replica-exchange molecular dynamics simulations for N1 and U14C, U14C[Formula: see text], U15A, and A17G mutants, ensuring extensive conformational sampling of both RNA and neomycin. The obtained neomycin association and binding paths, along with multidimensional free-energy profiles, revealed a two-step binding mechanism, consisting of conformational selection and induced fit. Neomycin binds to a preformed N1 conformation upon identifying a stable upper stem and U-turn motif in the riboswitch hairpin. However, the positioning of neomycin in the binding site occurs at different RNA-neomycin distances for each mutant, which may explain their different regulatory activities. The subsequent induced fit arises from the interactions of the neomycin's N3 amino group with RNA, causing the G9 backbone to rearrange. In the A17G mutant, the critical C6-A17/G17 stacking forms at a closer RNA-neomycin distance compared to N1. These findings together with estimated binding free energies coincide with experiments and elucidate why the A17G mutation decreases and U15A enhances N1 activity in response to neomycin.

Iron uptake pathway of Escherichia coli as an entry route for peptide nucleic acids conjugated with a siderophore mimic.

Bacteria secrete various iron-chelators (siderophores), which scavenge Fe3+ from the environment, bind it with high affinity, and retrieve it inside the cell. After the Fe3+ uptake, bacteria extract the soluble iron(II) from the siderophore. Ferric siderophores are transported inside the cell via the TonB-dependent receptor system. Importantly, siderophore uptake paths have been also used by sideromycins, natural antibiotics. Our goal is to hijack the transport system for hydroxamate-type siderophores to deliver peptide nucleic acid oligomers into Escherichia coli cells. As siderophore mimics we designed and synthesized linear and cyclic Nδ-acetyl-Nδ-hydroxy-l-ornithine based peptides. Using circular dichroism spectroscopy, we found that iron(III) is coordinated by the linear trimer with hydroxamate groups but not by the cyclic peptide. The internal flexibility of the linear siderophore oxygen atoms and their interactions with Fe3+ were confirmed by all-atom molecular dynamics simulations. Using flow cytometry we found that the designed hydroxamate trimer transports PNA oligomers inside the E. coli cells. Growth recovery assays on various E. coli mutants suggest the pathway of this transport through the FhuE outer-membrane receptor, which is responsible for the uptake of the natural iron chelator, ferric-coprogen. This pathway also involves the FhuD periplasmic binding protein. Docking of the siderophores to the FhuE and FhuD receptor structures showed that binding of the hydroxamate trimer is energetically favorable corroborating the experimentally suggested uptake path. Therefore, this siderophore mimic, as well as its conjugate with PNA, is most probably internalized through the hydroxamate pathway.

Structural dynamics influences the antibacterial activity of a cell-penetrating peptide (KFF)3K.

Given the widespread demand for novel antibacterial agents, we modified a cell-penetrating peptide (KFF)3K to transform it into an antibacterial peptide. Namely, we inserted a hydrocarbon staple into the (KFF)3K sequence to induce and stabilize its membrane-active secondary structure. The staples were introduced at two positions, (KFF)3K[5-9] and (KFF)3K[2-6], to retain the initial amphipathic character of the unstapled peptide. The stapled analogues are protease resistant contrary to (KFF)3K; 90% of the stapled (KFF)3K[5-9] peptide remained undigested after incubation in chymotrypsin solution. The stapled peptides showed antibacterial activity (with minimal inhibitory concentrations in the range of 2-16 µM) against various Gram-positive and Gram-negative strains, contrary to unmodified (KFF)3K, which had no antibacterial effect against any strain at concentrations up to 32 µM. Also, both stapled peptides adopted an α-helical structure in the buffer and micellar environment, contrary to a mostly undefined structure of the unstapled (KFF)3K in the buffer. We found that the antibacterial activity of (KFF)3K analogues is related to their disruptive effect on cell membranes and we showed that by stapling this cell-penetrating peptide, we can induce its antibacterial character.

Peptide nucleic acid conjugates and their antimicrobial applications-a mini-review.

Peptide nucleic acid (PNA) is a nucleic acid mimic with high specificity and binding affinity to natural DNA or RNA, as well as resistance to enzymatic degradation. PNA sequences can be designed to selectively silence gene expression, which makes PNA a promising tool for antimicrobial applications. However, the poor membrane permeability of PNA remains the main limiting factor for its applications in cells. To overcome this obstacle, PNA conjugates with different molecules have been developed. This mini-review focuses on covalently linked conjugates of PNA with cell-penetrating peptides, aminosugars, aminoglycoside antibiotics, and non-peptidic molecules that were tested, primarily as PNA carriers, in antibacterial and antiviral applications. The chemistries of the conjugation and the applied linkers are also discussed.

Conjugates of Aminoglycosides with Stapled Peptides as a Way to Target Antibiotic-Resistant Bacteria.

The misuse and overuse of antibiotics led to the development of bacterial resistance to existing aminoglycoside (AMG) antibiotics and limited their use. Consequently, there is a growing need to develop effective antimicrobials against multidrug-resistant bacteria. To target resistant strains, we propose to combine 2-deoxystreptamine AMGs, neomycin (NEO) and amikacin (AMK), with a membrane-active antimicrobial peptide anoplin and its hydrocarbon stapled derivative. The AMG-peptide hybrids were conjugated using the click chemistry reaction in solution to obtain a non-cleavable triazole linker and by disulfide bridge formation on the resin to obtain a linker cleavable in the bacterial cytoplasm. Homo-dimers connected via disulfide bridges between the N-terminus thiol analogues of anoplin and hydrocarbon stapled anoplin were also synthesized. These hybrid compounds show a notable increase in antibacterial and bactericidal activity, as compared to the unconjugated ones or their combinations, against Gram-positive and Gram-negative bacteria, especially for the strains resistant to AMK or NEO. The conjugates and disulfide peptide dimers exhibit low hemolytic activity on sheep red blood erythrocytes.

Varying molecular interactions explain aspects of crowder-dependent enzyme function of a viral protease.

Biochemical processes in cells, including enzyme-catalyzed reactions, occur in crowded conditions with various background macromolecules occupying up to 40% of cytoplasm's volume. Viral enzymes in the host cell also encounter such crowded conditions as they often function at the endoplasmic reticulum membranes. We focus on an enzyme encoded by the hepatitis C virus, the NS3/4A protease, which is crucial for viral replication. We have previously found experimentally that synthetic crowders, polyethylene glycol (PEG) and branched polysucrose (Ficoll), differently affect the kinetic parameters of peptide hydrolysis catalyzed by NS3/4A. To gain understanding of the reasons for such behavior, we perform atomistic molecular dynamics simulations of NS3/4A in the presence of either PEG or Ficoll crowders and with and without the peptide substrates. We find that both crowder types make nanosecond long contacts with the protease and slow down its diffusion. However, they also affect the enzyme structural dynamics; crowders induce functionally relevant helical structures in the disordered parts of the protease cofactor, NS4A, with the PEG effect being more pronounced. Overall, PEG interactions with NS3/4A are slightly stronger but Ficoll forms more hydrogen bonds with NS3. The crowders also interact with substrates; we find that the substrate diffusion is reduced much more in the presence of PEG than Ficoll. However, contrary to NS3, the substrate interacts more strongly with Ficoll than with PEG crowders, with the substrate diffusion being similar to crowder diffusion. Importantly, crowders also affect the substrate-enzyme interactions. We observe that both PEG and Ficoll enhance the presence of substrates near the active site, especially near catalytic H57 but Ficoll crowders increase substrate binding more than PEG molecules.

Trimeric Architecture Ensures the Stability and Biological Activity of the Calf Purine Nucleoside Phosphorylase: In Silico and In Vitro Studies of Monomeric and Trimeric Forms of the Enzyme.

Mammalian purine nucleoside phosphorylase (PNP) is biologically active as a homotrimer, in which each monomer catalyzes a reaction independently of the others. To answer the question of why the native PNP forms a trimeric structure, we constructed, in silico and in vitro, the monomeric form of the enzyme. Molecular dynamics simulations showed different geometries of the active site in the non-mutated trimeric and monomeric PNP forms, which suggested that the active site in the isolated monomer could be non-functional. To confirm this hypothesis, six amino acids located at the interface of the subunits were selected and mutated to alanines to disrupt the trimer and obtain a monomer (6Ala PNP). The effects of these mutations on the enzyme structure, stability, conformational dynamics, and activity were examined. The solution experiments confirmed that the 6Ala PNP mutant occurs mainly as a monomer, with a secondary structure almost identical to the wild type, WT PNP, and importantly, it shows no enzymatic activity. Simulations confirmed that, although the secondary structure of the 6Ala monomer is similar to the WT PNP, the positions of the amino acids building the 6Ala PNP active site significantly differ. These data suggest that a trimeric structure is necessary to stabilize the geometry of the active site of this enzyme.

Correction: Extracellular loops of BtuB facilitate transport of vitamin B12 through the outer membrane of E. coli.

[This corrects the article DOI: 10.1371/journal.pcbi.1008024.].

Stapled Anoplin as an Antibacterial Agent.

Anoplin is a linear 10-amino acid amphipathic peptide (Gly-Leu-Leu-Lys-Arg-Ile-Lys-Thr-Leu-Leu-NH2 ) derived from the venom sac of the solitary wasp. It has broad antimicrobial activity, including an antibacterial one. However, the inhibition of bacterial growth requires several dozen micromolar concentrations of this peptide. Anoplin is positively charged and directly interacts with anionic biological membranes forming an α-helix that disrupts the lipid bilayer. To improve the bactericidal properties of anoplin by stabilizing its helical structure, we designed and synthesized its analogs with hydrocarbon staples. The staple was introduced at two locations resulting in different charges and amphipathicity of the analogs. Circular dichroism studies showed that all modified anoplins adopted an α-helical conformation, both in the buffer and in the presence of membrane mimics. As the helicity of the stapled anoplins increased, their stability in trypsin solution improved. Using the propidium iodide uptake assay in Escherichia coli and Staphylococcus aureus, we confirmed the bacterial membrane disruption by the stapled anoplins. Next, we tested the antimicrobial activity of peptides on a range of Gram-negative and Gram-positive bacteria. Finally, we evaluated peptide hemolytic activity on sheep erythrocytes and cytotoxicity on human embryonic kidney 293 cells. All analogs showed higher antimicrobial activity than unmodified anoplin. Depending on the position of the staple, the peptides were more effective either against Gram-negative or Gram-positive bacteria. Anoplin[5-9], with a lower positive charge and increased hydrophobicity, had higher activity against Gram-positive bacteria but also showed hemolytic and destructive effects on eukaryotic cells. Contrary, anoplin[2-6] with a similar charge and amphipathicity as natural anoplin effectively killed Gram-negative bacteria, also pathogenic drug-resistant strains, without being hemolytic and toxic to eukaryotic cells. Our results showed that anoplin charge, amphipathicity, and location of hydrophobic residues affect the peptide destructive activity on the cell wall, and thus, its antibacterial activity. This means that by manipulating the charge and position of the staple in the sequence, one can manipulate the antimicrobial activity.

Insight into the Binding and Hydrolytic Preferences of hNudt16 Based on Nucleotide Diphosphate Substrates.

Nudt16 is a member of the NUDIX family of hydrolases that show specificity towards substrates consisting of a nucleoside diphosphate linked to another moiety X. Several substrates for hNudt16 and various possible biological functions have been reported. However, some of these reports contradict each other and studies comparing the substrate specificity of the hNudt16 protein are limited. Therefore, we quantitatively compared the affinity of hNudt16 towards a set of previously published substrates, as well as identified novel potential substrates. Here, we show that hNudt16 has the highest affinity towards IDP and GppG, with Kd below 100 nM. Other tested ligands exhibited a weaker affinity of several orders of magnitude. Among the investigated compounds, only IDP, GppG, m7GppG, AppA, dpCoA, and NADH were hydrolyzed by hNudt16 with a strong substrate preference for inosine or guanosine containing compounds. A new identified substrate for hNudt16, GppG, which binds the enzyme with an affinity comparable to that of IDP, suggests another potential regulatory role of this protein. Molecular docking of hNudt16-ligand binding inside the hNudt16 pocket revealed two binding modes for representative substrates. Nucleobase stabilization by Π stacking interactions with His24 has been associated with strong binding of hNudt16 substrates.

Vitamin B12 - Peptide Nucleic Acid Conjugates.

Vitamin B12 (cobalamin, Cbl) is an essential nutrient for all mammals and some bacteria. From a chemical point of view, it is a highly functionalized molecule, which enables conjugation at multiple positions and attachment of various cargoes. Both mammalian and bacterial cells have developed a specific transport pathway for the uptake of vitamin B12, and as a consequence, cobalamin is an attractive candidate for the delivery of biologically relevant molecules into cells. Indeed, hybrid molecules containing vitamin B12 in their structure have found various applications in medicinal chemistry, diagnostics, and biological sciences.Herein, we describe synthetic strategies toward the synthesis of vitamin B12 conjugates with peptide nucleic acid (PNA ) oligomers. Such short-modified oligonucleotides targeted at bacterial DNA or RNA can act as antibacterial agents if efficiently taken up by bacterial cells. The uptake of such oligonucleotides is hindered by the bacterial cell envelope, but vitamin B12 was found to efficiently deliver antisense PNA into Escherichia coli and Salmonella Typhimurium cells. This paves the way to the use of vitamin B12-PNA conjugates in antibacterial and diagnostic applications.Vitamin B12-PNA conjugates can be prepared via copper(I)-catalyzed azide-alkyne cycloaddition (CuAAC) that gives access to covalently linked hybrids or via connecting both building blocks by reduction-sensitive disulfide bridge. Both approaches require prior modification of vitamin B12 by incorporation of the azide moiety or via transformation of the native functional group into a moiety reactive toward thiols. Conjugation of vitamin B12 with PNA-tagged substrates efficiently furnishes designed conjugates.

Crowding affects structural dynamics and contributes to membrane association of the NS3/4A complex.

Using molecular dynamics simulations, we describe how crowded environments affect the internal dynamics and diffusion of the hepatitis C virus proteases NS3/4A. This protease plays a key role in viral replication and is successfully used as a target for antiviral treatment. The NS3 enzyme requires a peptide cofactor, called NS4A, with its central part interacting with the NS3 ß-sheet, and flexible, protruding terminal tails that are unstructured in water solution. The simulations describe the enzyme and water molecules at atomistic resolution, whereas crowders are modeled via either all-atom or coarse-grained models to emphasize different aspects of crowding. Crowders reflect the polyethylene glycol (PEG) molecules used in the experiments to mimic the crowded surrounding. A bead-shell model of folded coarse-grained PEG molecules considers mainly the excluded volume effect, whereas all-atom PEG models afford more protein-like crowder interactions. Circular dichroism spectroscopy experiments of the NS4A N-terminal tail show that a helical structure is formed in the presence of PEG crowders. The simulations suggest that crowding may assist in the formation of an NS4A helical fragment, positioned exactly where a transmembrane helix would fold upon the NS4A contact with the membrane. In addition, partially interactive PEGs help the NS4A N-tail to detach from the protease surface, thus enabling the process of helix insertion and potentially helping the virus establish a replication machinery needed to produce new viruses. Results point to an active role of crowding in assisting structural changes in disordered protein fragments that are necessary for their biological function.

Mutations of N1 Riboswitch Affect its Dynamics and Recognition by Neomycin Through Conformational Selection.

Short, structured fragments of non-coding mRNA may act as molecular switches upon binding specific ligands, regulating the translation of proteins encoded downstream this mRNA sequence. One switch, called riboswitch N1, is regulated by aminoglycosides such as neomycin. Nucleobase mutations in the apical loop, although distant from the binding pocket, significantly affect neomycin affinity and riboswitch regulatory efficiency. To explain this influence, we conducted molecular dynamics simulations using generalized replica exchange with solute tempering (gREST). Translation assay of a reporter protein in a yeast system shows that mutating A17 to G in the riboswitch apical loop reduces 6-fold the translation regulation efficiency of the mutant. Indeed, simulations of the unbound riboswitch show that G17 frequently stacks with base 7, while base 8 is stabilized towards the binding site in a way that it may interfere with the conformational selection mechanism and decrease riboswitch regulatory activity. In the riboswitch complexes, this single-point A to G mutation disrupts a strong hydrogen bond between nucleotides 5 and 17 and, instead, a new hydrogen bond between residue 17 and neomycin is created. This change forces neomycin to occupy a slightly shifted position in the binding pocket, which increases neomycin flexibility. Our simulations of the U14C mutation suggest that the riboswitch complex with neomycin is more stable if cytosine 14 is protonated. A hydrogen bond between the RNA phosphate and protonated cytosine appears as the stabilizing factor. Also, based on the cell-free translation assay and isothermal titration calorimetry experiments, mutations of nucleotides 14 and 15 affect only slightly the riboswitch ability to bind the ligand and its activity. Indeed, the simulation of the unbound U15A mutant suggests conformations preformed for ligand binding, which may explain slightly higher regulatory activity of this mutant. Overall, our results corroborate the in vivo and in vitro experiments on the N1 riboswitch-neomycin system, detail the relationship between nucleobase mutations and RNA dynamics, and reveal the conformations playing the major role in the conformational selection mechanism.

Vitamin B12-peptide nucleic acids use the BtuB receptor to pass through the Escherichia coli outer membrane.

Short modified oligonucleotides that bind in a sequence-specific way to messenger RNA essential for bacterial growth could be useful to fight bacterial infections. One such promising oligonucleotide is peptide nucleic acid (PNA), a synthetic DNA analog with a peptide-like backbone. However, the limitation precluding the use of oligonucleotides, including PNA, is that bacteria do not import them from the environment. We have shown that vitamin B12, which most bacteria need to take up for growth, delivers PNAs to Escherichia coli cells when covalently linked with PNAs. Vitamin B12 enters E. coli via a TonB-dependent transport system and is recognized by the outer-membrane vitamin B12-specific BtuB receptor. We engineered the E. coli ΔbtuB mutant and found that transport of the vitamin B12-PNA conjugate requires BtuB. Thus, the conjugate follows the same route through the outer membrane as taken by free vitamin B12. From enhanced sampling all-atom molecular dynamics simulations, we determined the mechanism of conjugate permeation through BtuB. BtuB is a ß-barrel occluded by its luminal domain. The potential of mean force shows that conjugate passage is unidirectional and its movement into the BtuB ß-barrel is energetically favorable upon luminal domain unfolding. Inside BtuB, PNA extends making its permeation mechanically feasible. BtuB extracellular loops are actively involved in transport through an induced-fit mechanism. We prove that the vitamin B12 transport system can be hijacked to enable PNA delivery to E. coli cells.

Conformational Changes of Anoplin, W-MreB1-9, and (KFF)3K Peptides near the Membranes.

Many peptides interact with biological membranes, but elucidating these interactions is challenging because cellular membranes are complex and peptides are structurally flexible. To contribute to understanding how the membrane-active peptides behave near the membranes, we investigated peptide structural changes in different lipid surroundings. We focused on two antimicrobial peptides, anoplin and W-MreB1-9, and one cell-penetrating peptide, (KFF)3K. Firstly, by using circular dichroism spectroscopy, we determined the secondary structures of these peptides when interacting with micelles, liposomes, E. coli lipopolysaccharides, and live E. coli bacteria. The peptides were disordered in the buffer, but anoplin and W-MreB1-9 displayed lipid-induced helicity. Yet, structural changes of the peptide depended on the composition and concentration of the membranes. Secondly, we quantified the destructive activity of peptides against liposomes by monitoring the release of a fluorescent dye (calcein) from the liposomes treated with peptides. We observed that only for anoplin and W-MreB1-9 calcein leakage from liposomes depended on the peptide concentration. Thirdly, bacterial growth inhibition assays showed that peptide conformational changes, evoked by the lipid environments, do not directly correlate with the antimicrobial activity of the peptides. However, understanding the relation between peptide structural properties, mechanisms of membrane disruption, and their biological activities can guide the design of membrane-active peptides.

Development of bis-ANS-based modified fluorescence titration assay for IFIT/RNA studies.

The interferon-induced proteins with tetratricopeptide repeats (IFITs) are a family of RNA-binding proteins that are very highly expressed during antiviral response of immune system. IFIT proteins recognize and tightly bind foreign RNA particles. These are primarily viral RNAs ended with triphosphate at the 5' or lacking methylation of the first cap-proximal nucleotide but also in vitro transcribed RNA synthesized in the laboratory. Recognition of RNA by IFIT proteins leads to the formation of stable RNA/IFIT complexes and translational shut off of non-self transcripts. Here, we present a fluorescent-based assay to study the interaction between RNA molecules and IFIT family proteins. We have particularly focused on two representatives of this family: IFIT1 and IFIT5. We found a probe that competitively with RNA binds the positively charged tunnel in these IFIT proteins. The use of this probe for IFIT titration allowed us to evaluate the differences in binding affinities of mRNAs with different variants of 5' ends.

Targeting Type II Toxin-Antitoxin Systems as Antibacterial Strategies.

The identification of novel targets for antimicrobial agents is crucial for combating infectious diseases caused by evolving bacterial pathogens. Components of bacterial toxin-antitoxin (TA) systems have been recognized as promising therapeutic targets. These widespread genetic modules are usually composed of two genes that encode a toxic protein targeting an essential cellular process and an antitoxin that counteracts the activity of the toxin. Uncontrolled toxin expression may elicit a bactericidal effect, so they may be considered "intracellular molecular bombs" that can lead to elimination of their host cells. Based on the molecular nature of antitoxins and their mode of interaction with toxins, TA systems have been classified into six groups. The most prevalent are type II TA systems. Due to their ubiquity among clinical isolates of pathogenic bacteria and the essential processes targeted, they are promising candidates for the development of novel antimicrobial strategies. In this review, we describe the distribution of type II TA systems in clinically relevant human pathogens, examine how these systems could be developed as the targets for novel antibacterials, and discuss possible undesirable effects of such therapeutic intervention, such as the induction of persister cells, biofilm formation and toxicity to eukaryotic cells.

Crowded environment affects the activity and inhibition of the NS3/4A protease.

Kinetic parameters characterizing the catalytic activities of enzymes are typically investigated in dilute solutions. However, in reality, these reactions occur in cells that, in addition to water and ions, are full of other macromolecules including proteins, nucleic acids, lipids, and metabolites. Such a crowded environment might affect enzyme-catalyzed reaction rates, so it is necessary to mimic the crowd in laboratory settings. We determined the effect of macromolecular crowders on the activity of the hepatitis C virus protease NS3/4A. As crowders we used polyethylene glycol (PEG), Ficoll, and bovine serum albumin. Using the fluorescence assay with a labeled peptide substrate, we found that the crowders affected the kinetics of the NS3/4A-catalyzed reaction differently. The Ficoll crowders increased and PEG decreased the initial and maximum reaction velocities. To explain the opposite effects exerted by PEG as compared to Ficoll, we performed molecular dynamics simulations of NS3/4A in explicit solvent and surrounded by its peptide substrates and PEG molecules. The simulations suggest both hydrophobic and polar/electrostatic interactions between PEG and NS3/4A with hydrogen bonds formed between PEG oxygens and NS3/4A amino acids rich in hydrogen bonds donors. The NS3/4A protease is a known target for telaprevir, an anti-viral drug. We found that Ficoll changes the inhibition constant for telaprevir suggesting that the effect of crowders should also be considered in inhibitor design.

Extracellular loops of BtuB facilitate transport of vitamin B12 through the outer membrane of E. coli.

Vitamin B12 (or cobalamin) is an enzymatic cofactor essential both for mammals and bacteria. However, cobalamin can be synthesized only by few microorganisms so most bacteria need to take it up from the environment through the TonB-dependent transport system. The first stage of cobalamin import to E. coli cells occurs through the outer-membrane receptor called BtuB. Vitamin B12 binds with high affinity to the extracellular side of the BtuB protein. BtuB forms a ß-barrel with inner luminal domain and extracellular loops. To mechanically allow for cobalamin passage, the luminal domain needs to partially unfold with the help of the inner-membrane TonB protein. However, the mechanism of cobalamin permeation is unknown. Using all-atom molecular dynamics, we simulated the transport of cobalamin through the BtuB receptor embedded in an asymmetric and heterogeneous E. coli outer-membrane. To enhance conformational sampling of the BtuB loops, we developed the Gaussian force-simulated annealing method (GF-SA) and coupled it with umbrella sampling. We found that cobalamin needs to rotate in order to permeate through BtuB. We showed that the mobility of BtuB extracellular loops is crucial for cobalamin binding and transport and resembles an induced-fit mechanism. Loop mobility depends not only on the position of cobalamin but also on the extension of luminal domain. We provided atomistic details of cobalamin transport through the BtuB receptor showing the essential role of the mobility of BtuB extracellular loops. A similar TonB-dependent transport system is used also by many other compounds, such as haem and siderophores, and importantly, can be hijacked by natural antibiotics. Our work could have implications for future delivery of antibiotics to bacteria using this transport system.